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Abstract 

We discuss leptogenesis within a TeV-scale inverse seesaw model for neutrino masses where 
the seesaw structure is guaranteed by an 50(10) symmetry. Contrary to the TeV-scale type-I 
gauged seesaw, the constraints imposed by successful leptogenesis in these models are rather 
weak and allow for the extra gauge bosons Wji and Z' to be in the LHC accessible range. The 
key differences in the inverse seesaw compared to the type I case are: (i) decay and inverse decay 
rates larger than the scatterings involving extra gauge bosons due to the large Yukawa couplings 
and (ii) the suppression of the washout due to very small lepton number breaking. 
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I. INTRODUCTION 

One of the attractive features of the seesaw mechanism for neutrino masses [l| 
is that it provides a way to understand the origin of matter in the Universe via leptogen- 
esis {2! (for a recent review, see Ref. jsj). In the vanilla framework of leptogenesis where 
right-handed (RH) neutrino masses are hierarchical, it is well known that the lightest RH 
neutrino needs to be rather heavy, around 10^ GeV or higher ^|. These scales are however 
beyond the reach of collider experiments, e.g. the CERN Large Hadron Collider (LHC). 
On the other hand, from the point of view of the seesaw model itself, one can envisage 
the new physics scale to be anywhere between TeV to 10^^ GeV. It is well known that 
the scale of leptogenesis can be lowered to the TeV scale if one allows the RH neutrinos 

n 

to be quasi-degenerate |5|. However, first the quasi-degeneracy should be motivated, and 
second if the scale of the RH neutrinos is to be explained by the breaking of some gauge 
symmetry, what is the impact of the latter on leptogenesis? 

Two classes of seesaw models are of interest in this connection: the usual 
type-I seesaw and the inverse seesaw 6|. In both classes of models, a higher gauge 
symmetry, e.g. B — L, is usually called for to make the model "natural". In addition to 
providing a compelling reason for the inclusion of the RH neutrinos to guarantee anomaly 
cancellation, in the type-I case it can be used to understand why the seesaw scale is so 
much lower than the Planck scale, whereas, in the inverse seesaw case, it stabilizes the 
zeros in the (z/, A^, S) mass matrix that leads to the doubly-suppressed seesaw formula. 

An attractive gauge symmetry that embeds the B — L symmetry and also 
provides a way to understand the origin of parity violation in low-energy weak interactions 
is the Left-Right (LR) gauge group S't/(2)ixS'f/(2)ijxt/(l)B_i:, Q. An important question 
that arises in these models is: What is the scale of parity invariance? In particular, if it is 
in the TeV range and if at the same time leptogenesis generates the desired matter-anti- 
matter asymmetry, then the LHC could be probing neutrino mass physics as well as shed 
light on one of the deepest mysteries of cosmology. 

Since Sakharov's out-of-equilibrium condition jsj] must be satisfied in order to 
generate a baryon asymmetry, the existence of new interactions inherent to the LR mod- 
els make it a nontrivial task to check whether a TeV-scale Wr is indeed compatible with 
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leptogenesis as an explanation of the origin of matter. Specifically, the efficiency of lep- 
togenesis crucially depends on the number of RH neutrinos that decay out of equilibrium 
to produce a leptonic asymmetry. This number is set by two things: First, it depends 
on the relative magnitudes of the decay rate and the (C7P-conserving) gauge scattering 
rates of the RH neutrino, since this can lead to a dilution of the number of "useful" RH 
neutrinos. Second, the washout processes, primarily inverse decays, should drop out of 
equilibrium early enough, otherwise the number of RH neutrinos gets suppressed at an 
exponential rate. 

These issues have been analyzed for the type-I case within LR symmetric 
models {q] as well a.s B — L models 10|. It was found that for the full LR models with 
TeV-scale parity restoration and RH neutrino masses, gauge scattering rates induced by 
Wr exchange largely dominate the decay and inverse decay rates because the Yukawa 
couplings are small for the standard type-I seesaw at the TeV scale. These facts lead 
to a huge dilution of the number of RH neutrinos which decay out of equilibrium and 
in a CP asymmetric manner. Moreover, the gauge scattering interactions also wash out 
lepton number at a very large rate, much larger than the inverse decays. Altogether, 
these two effects lead to a very stringent constraint on the mass scale of Wfi for successful 
leptogenesis, Mw^^ > 18 TeV |9|, which would imply that the discovery of a Wr at the LHC 
is incompatible with thermal leptogenesis as the origin of matter. On the other hand, in 
the case of a simple B — L theory, successful leptogenesis only implies that Mz' > 2.5 TeV 
in the "collider-friendly" region of parameter space where the RH neutrino mass is less 



than half the Z' mass 
energy observations 



10 



p. We note that there exist bounds on the Wr mass from low 

12| and they allow Wr mass to be as low as 2.5 TeV. 

In this paper, we have analyzed the leptogenesis constraints on the receritly 

proposed TeV-scale LR model within a unified supersymmetric 5*0(10) framework 13| 

where neutrino masses arise from an inverse seesaw mechanis Two features distinguish 

the inverse seesaw mechanism from the type-I seesaw: (i) the Dirac Yukawa couplings of 

the RH neutrino can be much larger (~ 10^^ — 10"^) than for the type-I case (where 

they are typically of order ~ 10^^ for TeV-scale RH neutrino masses) and (ii) the lepton- 

^ For a discussion of low scale leptogenesis in an S0{1Q) model where only the doubly charged Higgs 

boson is in the TeV range, see Ref. 11 1. 
^ For other low-scale leptogenesis scenarios in inverse-seesaw-related frameworks, see Ref. 

3 
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number- violating parameter (the Majorana mass fi of the left-right singlet lepton S, which 
measures the "pseudo-Diracness" of N) is much smaller than the Dirac mass of A^. As 
a result, first, the decay rate of can be much larger than the Wr exchange scattering 
rate at the baryogenesis epoch, and second, the wash-out processes are suppressed by the 
small Majorana mass fi. Consequently, we find that both the Wji and Z' can be in the 
TeV range and hence accessible at the LHC This is the main result of our paper, and it 



should make the case for searching the Wn and Z' at LHC stronger |15 |. 

This paper is organized as follows: in Section II, we summarize the LR inverse 
seesaw model and give the Dirac Yukawa couplings as well as the various lepton-number- 

n 

violating parameters as constrained by SO (10) unification p^. In Section III, we present 
a generic discussion of leptogenesis in this class of inverse seesaw models; in Section IV, 
we present the numerical results for our model. Finally, we summarize our findings and 
conclude in Section V. In Appendix A, we present a new scenario for gauge coupling 
unification (different from that discussed in Ref. jisl]) in these models where the relative 
magnitudes of Wr and Z' masses can be unrelated. In Appendix B, we give the analytical 
expressions of the CP-asjTiimetry in the inverse seesaw model for some special cases. 

II. LEFT-RIGHT INVERSE SEESAW PARAMETERS IN 50(10) 



The implementation of the inverse seesaw mechanism p] requires, in addition 
to the usual Standard Model (SM) singlet RH neutrinos A^, {i = 1,2,3 for three genera- 
tions) as in the typical type-I seesaw, three extra SM gauge singlet fermions Si coupled 
to the RH neutrinos through the lepton-number-conserving couplings of the type NS, 
while the traditional RH neutrino Majorana mass term is forbidden by the lepton number 
symmetr}0. In the low energy theory, dominant lepton number breaking arises only from 
the self-coupling term SS. The neutrino mass Lagrangian in the flavor basis is given by 

C^^s = u^C-^MdN + N'^C-^M^S + i^^C-V^ + h-c, (1) 

where /i is a complex symmetric 3x3 mass matrix containing all the lepton-number- 
violating parameters, and M£, and M^v are 3x3 mass matrices representing the Dirac 
^ If we include higher dimensional terms in the theory, they can induce an NN Majorana mass term but 
its magnitude is of order w^^/Afpi ^ 10^^"^ GeV and is too small to affect our discussion. 
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mass terms in the v-N and N-S sectors, respectively. In the basis {z/, A^, S*}, the full 9x9 
neutrino mass matrix is then given by 

^0 Md ^ 
M^= Ml Mn ■ (2) 

The lepton-number-violating entries in the /i matrix have to be much smaller than the 
Dirac neutrino masses in order to fit the light neutrino masses, as observed in neutrino 
oscillation experiments. In fact, the light neutrino mass matrix can be cast in a seesaw-like 
form in the limit fi <^ Md <^ M^: 

~ MoM^^fi (m^) Mj; = F^F^, (3) 

to leading order in F = M^M^^. As expected, in the limit /i — 0, which corresponds to 
unbroken lepton number, we recover the massless neutrinos of the SM. We note that this 
smallness of the yU-parameter peculiar to the inverse seesaw models allows for a neutrino 
mass fit even with TeV-scale RH neutrino mass and large Dirac mass terms. Theoretically, 
smallness of the u-term could be explained in extra dimensional brane world models if 

n 

the lepton number is broken in a separate brane from the standard model brane 1161 . 

n 

As shown in Ref . , in order to embed a TeV-scale inverse seesaw mechanism 
into a generic SO{10) model, we need to break the B — L gauge symmetry by 16-Higgs 
fields at the TeV scale, whereas the S'O(IO) symmetry is broken down to the LR symmetric 
gauge group 5f/(3)c x SU{2)l x SU{2)r x U{1)b-l at the GUT-scale by 45 and 54-Higgs 
fields; finally, the SM symmetry is broken at the weak scale by 10-Higgs fields. As in 
the usual 5*0(10) models, the three generations of quark and lepton fields are assigned to 
three 16-dimensional spinor representations, and correspondingly, we add three 5*0(10) 
singlet matter fields Ij (they can be identified with the Si fields above) to implement the 
inverse seesaw mechanism. 

As discussed in Ref. [l^, we need at least two 10// fields to have a realistic 
fermion mass spectrum; we also need two 45// fields, one for symmetry breaking at the 
GUT scale and another to give rise to the vectorlike color triplets at the TeV-scale as 
required by coupling unification constraints. Similarly, we need only the SU{2)r doublet 
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fields of 16h and no SU (2) l fields for unificatioij^- With this minimal set of Higgs fields, 
the most general Yukawa superpotential is given by 

Wy = KijWdGAOH^ + ^16a6,10H„45H45^ + ^ 16,16, 10^,^45^, 

+/^,16a,T6H,+/i.,ia„ (4) 

where the first term is the usual Yukawa coupling term, the second and third terms are 
higher-dimensional terms, and the last two terms give rise to the inverse seesaw mecha- 
nism. As already pointed out in Ref. 13|, it is sufficient to keep only one of the higher- 
dimensional operators, usually the 16 ■ 16 ■ 10 ■ 45 ■ 45' term, whose fully antisymmetric 
combination acts as an effective 126^;^ operator, in order to obtain a realistic fermion mass 
spectrum at the GUT scale, and hence for simplicity, we will assume all the /'-couplings 
to be zero; keeping this term does not affect our discussion below [fj. 

The B — L symmetry is broken when the 16//-field acquires a vacuum expec- 
tation value (VEV) and the N-S sector RH neutrino mass matrix is given by 

Mn,, = ^fi/^, (5) 

where vr is the VEV of 16h and is of order TeV for the low-scale B — L breaking models 
considered here. All the other fermion masses are generated when the SM symmetry is 
broken at the weak scale by the 10// VEVs. We consider here only the model (A) of 



Ref. 



13| where the VEV patterns of the two 10// fields are given by 

^i)= \ °V {^2)={'' "I, (6) 






and the fermion mass matrices are given by 

Mu = hu + /, Md = hd + /, 

Me = /i,-3/, MD = k-3f, (7) 

* An alternative choice of Higgs fields which also consistently leads to coupling unification in this scenario 

is presented in Appendix A. 
^ The /' term has two effective contributions - one of 10 -Higgs type and another of 120-Higgs type. The 

effective 10 coupling can be absorbed into the first term, and since the 120 coupling is antisymmetric 

in generation indices, it only contributes to the off-diagonal elements in fermion mass matrices. Hence, 

a non-zero /' coupling could only slightly modify the specific structure of the Dirac neutrino mass 

matrix, without changing any of the main results of the paper. 
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where in the notation of Ref. hu,d = Vu,dhu,d and / = Vufu = Vdfd- Using the 
renormahzation group evolution of the fermion masses in the LR model, we obtain the 
GUT-scale fermion masses starting from the experimentally known weak scale values, and 
using these mass eigenvalues, we obtain a fit for the Yukawa coupling matrices at the GUT 
scale, from which we can get the structure of the Dirac neutrino mass matrix. Here, as 
an example, we quote the result for tan/3 = Vu/vd = 10 jl3 |: 

^ 0.0111 0.0384- 0.01031 0.038 - 0.44331 ^ 
Md= 0.0384 + 0.01031 0.2928 1.8623 + 0.00021 GeV. (8) 

^ 0.038 + 0.44331 1.8623 - 0.00021 77.7573 y 
With this Dirac neutrino mass, we can easily fit the observed neutrino oscillation data 
by fixing the singlet mass matrix /i in Eq. ([3]). As an example, for a normal hierarchy of 
neutrino masses, and assuming a diagonal structure for the RH neutrino mass matrix Mm 
with eigenvalues (3.5, 3, 1) TeV, we can fit the observed 2a neutrino oscillation data [l7| 
for the following choice of /i: 

^-17.4226 + 0.30981 2.1033 - 0.05901 -0.0136 + 0.02881 ^ 
^^= 2.1033 - 0.05901 -0.2585 + 0.00971 0.0016 - 0.00351 GeV. (9) 

^-0.0136 + 0.02881 0.0016 - 0.00351 3.6 x 10^^ + 4.6 x lO'^i y 



III. LEPTOGENESIS IN LEFT-RIGHT INVERSE SEESAW MODELS 



In this section we summarize the main features of leptogenesis within the class 
of LR inverse seesaw models discussed above. We also wish to note that while we have 
used the S'0( 10) framework to make the results definite and somewhat more predictive, 
our discussion applies also to the case with TeV-scale Left-Right symmetry without grand 
unification. In what follows, we will partially follow the discussion presented in Ref. |l^ . 

As discussed in the introduction, a crucial difference of the inverse seesaw 
from the usual seesaw is the dependence on a new mass matrix, /i, which can lead to the 
result that no matter what the ratio of the mass scales Md/M^ is, the lightness of the 
left-handed neutrinos can always be explained by small /i entries. In other words, the 
inverse seesaw makes it possible to have at the same time large Dirac masses and low, say 
TeV-scale, RH neutrino masses, and it still can explain why neutrinos are light. This is 



directly connected to the fact that, in the hmit — )■ 0, lepton number is conserved, and 
therefore neutrino masses vanish, as in the SM. This is a crucial difference from the case 
of TeV scale type I seesaw. 

An interesting question to ask is how leptogenesis is affected by this distinctive 
feature of inverse seesaw models. We expect that the lepton-number-violating washout 
will go to zero in the limit of vanishing fi. As a matter of fact, as explicitly shown in 



Ref. 



18|, the all-important AL = 2 washout process — > vanishes as 6f, with 



im^^f^ (10) 

where Fj is the total decay rate of Ni into lepton and Higgs (and antiparticles), and 
Mij are the masses of the quasi-Dirac RH neutrino pair Nij (with Fj ~ F^). Note that 
we denote by Mj {i = 1, . . . , 6) the heavy neutrino mass eigenvalues. As shown in the 
Appendix, the leading order contribution to the mass splitting for each quasi-Dirac pair 
comes from the diagonal elements of the /i matrix. Therefore, as expected, the washout 
tends to zero in the limit of vanishing /x. The suppression of the washout can be shown to 
occur through the destructive interference of one member of a quasi-Dirac pair with the 
other 181]. It is instructive to show numerically how the washout is kept under control in 



this family of models with more than one pair of RH neutrinos. The washout parameter 
Ki is defined as 

_ {hh%vl 

' H{z = l) m.M, ' 
where H(z) is the usual Hubble expansion rate: H{z) ^ 1.66y^^2]vf^, ^* — ^-^^ ^ 
10~^ eV|j. Plugging in numbers, we find that with Yukawa couplings of order 10~^ and 
a RH neutrino mass of order 1 TeV, the washout parameter K is of order 10^^, which is 
huge! However, the suppression of the washout is also very large, being proportional to 
with 5^1 due to the smallness of fi, as required to get the right scale for the light 
neutrinos. Specifically, for the example of Eq. we find that 6 ~ 10^^ and therefore 
the effective washout parameter K'^^ ~ 6^K ~ 100, which is reasonably small. 
^ Note here that for simphcity, we have assumed the SUSY breaking scale to be above the hghtest 

RH neutrino mass so that only the SM degrees of freedom are in relativistic thermal equilibrium, i.e. 

g^, ~ 106.75. However, the main results of this paper remain unchanged irrespective of the sparticle 

spectrum chosen. 
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In the LR model we are considering, there are other processes contributing to 
the washout of lepton number, for instance, NjieR ^ u^d^. More precisely, this process 
destroys RH lepton number, but in the temperature range of interest to us (TeV scale) 
every individual RH lepton flavor equilibrates with the LH lepton flavor one, thanks to 
the Yukawa interactions. Does this process also turn off in the limit of lepton number 
conservation? It can be easily shown that, including the production of the RH neutrino 
by an inverse decay, followed by the scattering process mentioned above, there is also a 
destructive interference within the quasi-Dirac pair which leads exactly to the same kind 
of 5^-suppression as for the process — > I^'^ . 

Another feature of inverse seesaw models is that they typically lead to lepton 
flavor equilibration [2^ because of the large Yukawa couplings. More precisely, it can be 
shown that the process £q,<I> ^ which does not change lepton number, but changes 
lepton flavor, is deep in thermal equilibrium for the TeV temperatures (see, for instance, 
Ref. 2l(]). Consequently, the Boltzmann equations for leptogenesis can be written as only 
one equation for the sum of the lepton flavors 20]. In other words, flavor effects |3| are 
not important in our framework. 

Putting together all the qualitatively important effects discussed above and 
solving the relevant set of Boltzmann equations, one can derive the following expression 
for the efficiency factor (see, for instance, Ref. 22|): 

,rfiV^^(/) D{K„z') 



dz' D{K„ z') + DwA^') + ^Sz'N^iz') + SwA^') 



xexp 



(12) 



where z = Mi/T, Nx is the number density of X over the relativistic number density 
of RH neutrinos, and D, S, W denote the various decay, scattering and washout terms 
respectively, defined in Section IV. Note that the expression above assumes dNj^Jdz ~ 
dN^^Jdz with A''^ = ^^^^^ii^z), where lC2{z) is the modified Bessel function of the 2nd 
kind. This is a very good approximation in our model (with large Yukawa couplings). 
Note also that we are neglecting spectator processes 231] and AL = 1 scatterings involving 
the Higgs, which are both expected to lead to order one corrections. 
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The final baryon asymmetry can be conveniently written as 

r]B ^ 10'^^ eiKi{z oo) , (13) 

i 

where the dilution factor 10~^ takes into account the fraction oi B — L asymmetry con- 
verted into baryon asymmetry by sphaleron processes and also the dilution due to photon 
production from the onset of leptogenesis till recombination, is the CP-asymmetry 
generated by the decay of Ni into any lepton flavor and is given by 2^ 

where f" is the L-violating self-energy and vertex loop factoiQ. In the quasi-degenerate 
limit of the (^, j) pair, we have 

M2 - Mf 

Note that Eq. f |T^ was derived assuming heavy neutrino mass eigenstates. 
Therefore, it is necessary to make a basis transformation from the "flavor" basis where 

/ 



M 



RH 



(16) 



Mn 

to the diagonal mass basis with real and positive eigenvalues Mj {i = 1, 2, . . . , 6), grouped 
into three quasi-degenerate pairs with mass splittings in each pair of order /i^fc {k = 1, 2, 3). 
Analytically, the exact diagonalization of the full 6x6 mass matrix Ainn to get a closed 
form expression for the Yukawa couplings, hia, in terms of the known parameters, namely 
M/), Mjv and /i, is extremely involved. In Appendix B, we show the analytical expressions 
up to first order in fi for some simpler cases with only two quasi-Dirac pairs and show 
explicitly that the CP-asymmetry indeed vanishes in the L-conserving limit /i — )■ 0, as 
expected. For the general case with three quasi-Dirac pairs, we numerically evaluate the 
CP-asymmetry in the next section. We note that the three-pair case reduces to the two- 
pair case discussed in Appendix B if one of the masses is much heavier than the other 
two and hence decouples from the rest. 



Note that the L-conserving self-energy contribution vanishes when one sums over flavor. 
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IV. RESULTS 



As noted in Section II, the Yukawa couplings are fixed by the 5*0(10) symme- 
try. The fi matrix can be deduced from the knowledge of the light neutrino masses and 
mixing angles as a function of the RH neutrino mass matrix M^v, which can be taken to 
be diagonal without loss of generality. Varying the RH neutrino mass eigenvalues input 
then leads to different /i matrices, keeping the light neutrino mass matrix given by Eq. (|3]) 
such that its mass eigenvalues and mixing angles are within 2a of the observed values. 
Once we know the explicit form of the RH neutrino mass matrix given by Eq. (fT6ll . we 
can define the quasi-Dirac pairs by transforming to a basis in which this mass matrix is 
diagonal with real and positive eigenvalues. We then calculate the CP-asymmetry and 
efficiency factors for the decay of the lightest RH neutrino pair and scan the parameter 
space to match the calculated baryon asymmetry (using Eq. f fT3|) ) with the observed 68% 



CL. value, tjb = (6.2 ± 0.15) x 10^^" 25|]. Note that we only consider the asymmetry 
generated by decay of the lightest RH neutrino pair as the asymmetry generated by the 
heavy pairs is washed out very rapidly (due to large exponential suppression), and for 
these washouts not to affect the asymmetry generated by the lightest pair, we require the 



lightest pair to be at least 3 times smaller than the next heavy pair 



26|. 



To calculate the efficiency factor given by Eq. ( [T2|) . we first write down the 



22|: 



thermally averaged rates for N ^ decay and the corresponding inverse decay 

WiY,{K,,z) = h<Xi{z)z\ (17) 

with /Cq, a = 1,2 denoting the modified Bessel function of the zth type. The thermally 
averaged rate Dw^ for the V^R-mediated A^-decay, — )• iqnq'ji, is given by 

D„Az) = ^ (18) 

where rf^ is the RH neutrino equilibrium number density, rf^{z) = ^n^{z)N'^^{z) with 
n-y = ^^^T^, and 7Ar is the reaction density: 

tIT"' = "5J(^)^rr«', (19) 

11 



where rj^^'' is the total three body decay width of A^, given by 



^ N 



^9% 



ds- 



(20) 



with the Wr total decay width = ^Mwji- 

The various scattering rates Swji,z' appearing in Eq. (fT2|) are also defined as 
in Eq. ( fTSl) where the corresponding scattering reaction density is related to the reduced 



cross section as follows (see, for instance, Ref. 



27|): 



■y{ab ^ cd) 



(21) 



with X = s/Mfj and the threshold value Xthr = T72-niax[(ma + m;,)^, {rric + "^d)^J- The 



reduced cross sections for various Wr exchange diagrams were computed in Ref. |9|]: 

a (r) - 94 (1 - + 2^^) .021 



AS-nx 



X 



Mir 
hp 



9gji ^^{x + u){x + u-l) 



iTTX Jl-x 



U 



Mir 



O'Nd 



R^ejiUR 



X] 



9g% M\ 



N 



X] 



Stt Ml 



Wr [X + 



Ml, 

^R 



(23) 
(24) 



Here we have ignored the t-channel process A^A^ — ^l as the rate for this process falls off 
very rapidly for the region of interest, viz. z > 1 jof. 

The reduced cross section for the Z' exchange diagram is given by 



X 



13g' 



B-L 



yjx{x — A)-- 



X 



Ml, 



(25) 



with the total Z' decay width 

Tz' = 



2Atx 



Mz' 



13 + 3 1 



4M|' 
Ml, 



3/2 



(26) 



Before calculating the efficiency factor, it is instructive to compare all the 
reaction rates appearing in Eq. (fT2l) to get a clear idea of various contributions. As 
an illustration, we consider the case with the RH Majorana neutrino mass eigenvalues 
(3.5,3,1) TeV (as in Eq. (|9])). The fiavor-summed washout parameter for the decay of 
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the lightest quasi-Dirac pair in this case is given by ^ 4 x 10 whereas the effective 
washout parameter is given by Kf^ ~ (^l-^s — 168 which is reasonable. For comparison, 



the corresponding values for the two heavy pairs are Kf^ ^ = 8 x 10^ which, when 
exponentiated in the washout term in Eq. (JT2l) . leads to a huge suppression, thus making 
the efficiency in those channels practically negligible. Hence, from now on, we will consider 
the decay of only the lightest pair. 




100 
0.01 



FIG. 1. Various decay and scattering rates (thermally averaged) as a function of z = M^^/T 
for a particular choice of RH neutrino masses, {M^-^^, AIj\f^, M^.^) = (3.5,3, 1) TeV and = 
Mz' = 2 TeV. The yellow shaded region is where the asymmetry is generated. 

In Fig. [H we show the various thermally averaged decay and scattering rates 
as a function of ^ = M^vg/T, for the above choice of the RH neutrino masses and for 
= Mz' = 2 TeV. The yellow shaded region shows the asymmetry production time. 
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We note 



approximately when zb — 2<z<zb + 2 with zb — 2 + 4^}{0.i3^-2.5/k ^ g 
that in this range, the — )■ decay rate, -D^$, dominates over the three-body decay rate 
as well as all the scattering rates by several orders of magnitude. Hence in the efficiency 
factor, Eq. ( !T2|) . the dilution term D/{D + S) is very close to unity and is essentially 
independent of Mw^^ and Mz'. The enhanced — > decay rate is due to the large 
Yukawa couplings in the inverse seesaw scenario. We also note that as the ly/j- mediated 
* Here we have assumed that the production of asymmetry stops immediately after the temperature drops 
below the sphaleron freeze-out temperature, Tgph — 130 GeV for a Higgs mass mu = 120 GeV y . 
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three-body decay rate is much smaller than the — > decay rate, the washout term 

in Eq. ( !T2|) arising due to the process — )• — )• Iqq' which is proportional to 

the branching ratio of N ^ iqq' will be suppressed compared to the inverse decay term 

WiD- Thus we find that the efficiency factor is also essentially independent of both Wr 

and Z' masses for a wide range of parameter space. Of course, the Wn and Z' scattering 

terms will start to dominate for very low values of their masses; however, we estimated 

this lower bound to be well below the current collider bounds on Mw„ and Mz> which are 
□ 

roughly a TeV or so [30|. 
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FIG. 2. Correlation between the baryon asymmetry and the lightest RH neutrino mass for 
various heavy mass pairs (M^Vi , M^^) in TeV. The yellow shaded region is the observed value of 
rjB within 2a C.L.. 



Fig. [2] shows the baryon asymmetry rjB as a function of the lightest RH neutrino 
mass M^Vg for different choices of the heavy RH neutrino masses (M^r^, Mjvj) from 1.5 — 3 
TeV. The calculated value of rjB is to be compared with the observed 68% C.L. value, 
Vb = (6.2 ± 0.15) X 10^^*^ 25|. It is clear from the ffgure that for a given heavy mass 
pair, there is a narrow range of values allowed for M^,^ satisfying the observed baryon 
asymmetry (the yellow shaded region). We note that for ffxed M^r^, the allowed range 
of Mjvg decreases with increasing M^^, while for ffxed Mjv^, the allowed range of Mjy^ 
increases with increasing M^r^. Also note that when the heavy pairs have degenerate 
mass, the baryon asymmetry gets suppressed (e.g. the lower two hues in Fig. [2]) due to 



14 



the suppression in the CP-asymmetry. Finally, we note that for a given set of heavy mass 
pairs, r]B oc M^^. 

Fig. [3] shows the correlation between the efficiency factor and the flavor- 
summed CP asymmetry for various channels. The different lines correspond to different 
values of the heavy mass pair, {M^-^^, Mj^i^) TeV, starting from (3,1.5) TeV at the top 
to (3,3) TeV at bottom. We note that for fixed M^r^, the lines move down as we increase 
MjVj while for fixed Mn2, they move up with increasing Mjvj. The yellow shaded region 
shows the observed value of tib which is essentially the product of k, and e, summed over 
all pairs. As we have pointed out earlier, only the lightest pair contribution is significant, 
while the efficiency is too small for the other two pairs. 




0.00002 0.00004 0.00006 0.00008 0.0001 



FIG. 3. Correlation between the efficiency factor and the flavor-summed CP-asymmetry for the 
lightest pair for various values of the heavy mass pair (Mtv^ , M^Tj ) in TeV. The yellow shaded 
region corresponds to the observed value of rjB within 2a C.L.. 

V. CONCLUSION 

In summary, we have shown that a TeV-scale Left-Right symmetry can be 
compatible with the understanding of the origin of matter via leptogenesis provided small 
neutrino masses are understood using the inverse seesaw mechanism. A crucial feature 
of this mechanism is that the magnitude of the lepton-number-breaking Majorana mass 
term is directly proportional to the neutrino mass, rather than inversely as in the usual 
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type-I seesaw framework. This allows the Yukawa couplings that generate the Dirac mass 
for the neutrinos close to one, even with TeV-scale RH neutrinos. These two facts help 
to keep the wash-out of the generated lepton asymmetry under control, and thus explain 
the origin of matter while keeping both the Z' and in the TeV range. The results of 
this paper should provide new motivation for searching for the and the Left-Right Z' 
at the LHC As has been already emphasized in literature 1^, the signal for the inverse 
seesaw with would be presence of trilepton final states with missing energy. 
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Appendix A: A New Coupling Unification Scenario 



It was shown in Ref. 13| that the 5*0(10) embedding of TeV-scale inverse 
seesaw mechanism discussed in this paper is consistent with gauge coupling unification. 
In particular, it was shown that for the symmetry breaking chain 

5-0(10) ^ 3,2i2^lB-L (SUSYLR) ^ 3,2ily (MSSM) ^ 3,2ay (SM) ^ 3,1q, 

and for TeV-scale SUSY and B — L breaking, it is possible to obtain gauge coupling 
unification with two IOh bidoublets, two RH 16h doublets and one 45h color triplet. In 
this section, we show that there is an alternative choice of Higgs fields which also leads 
to coupling unification with TeV-scale Wr and Z'. If we add one set of SU{2)r triplets, 
A(l, 1,3,0), coming from the 45// field, then we can also achieve unification with the 
same field content as above except that we need only one SU{2)r doublet. The SUSYLR 
/9-function in this model is given by 

/ 3 3 \ /23 \ 

hi = (^10 + -riL + -ur, nio + n^, mo + n^j + 2, -2j = (^y , 2, 3, -2j {i = I^-l, 2^, 2^, 3,) 
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for riio = 2, rii = and n/j = 1. As shown in Fig. |U we achieve unification at Mq — 
4 X 10^6 GeV with ~ 20.3. 

80 
70 
60 
50 

-^B 40 
30 
20 
10 


2 4 6 8 10 12 14 16 18 

t [logio(GeV)] 

FIG. 4. Gauge coupling unification in the new SO{10) scenario with one SU{2)f^ doublet and 
one SU{2)r triplet. 

This model has two nice features over the previous one: (i) all the Higgs fields 
required for unification are connected to breaking of separate gauge symmetries and there 
is no arbitrariness in the number of fields, and (ii) the presence of the SU{2)r triplet 
enables us to decouple the mass scales Mwji and Mz' which are otherwise related in usual 
Left-Right models with Mz' > My/^- 

Appendix B: Analytical Expression for the CP-Asymmetry 

The full Majorana mass matrix in the flavor basis (A^j, Si) is given by 

/ Mtv 

M = 

where for i = 1,2, 3, both M^r and /i are 3x3 symmetric matrices. The Yukawa Lagrangian 
in this basis is given by (with i,j = 1, 2, 3) 

Cy = ViaN.^n^ + M^^^NjC-^S, + ]^^i,,SlC-^S, + h.c. (B2) 
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In order to calculate the CP asymmetry in this framework, it is more convenient 
to work in the basis in which the RH Majorana neutrino mass matrix is diagonal with real 
and positive eigenvalues. The Lagrangian in this basis is given by (with i = 1, 2, • • • , 6) 

Ch = hiJlMa + ^MiNfC-^Ni + h.c. (B3) 

Analytically, the exact diagonalization of the full 6x6 mass matrix M. is extremely 
involved and we cannot obtain a closed form expression for the CP-asymmetry in this 
case. However, we can study the dependence of the small L-violating parameter in 
some special cases, viz. when the //-matrix is completely diagonal or completely off- 
diagonal, as in these cases the Majorana mass matrix reduces to a block diagonal form. 
In this section, we derive the analytical expression for the CP- asymmetry in these two 
limits and for two sets of RH neutrinos, i.e. for {Ni,Si) with i = 1,2. The i = 3 case 
reduces to this limit if one of the masses is much heavier and hence decouples from the 
other two. 

We consider the 4x4 version of the Majorana mass matrix M.: 

KA ( ° ^^-^ 

M4x4 = 

y ^N2x2 A*2x2 

where without loss of generality we choose the mass matrix to be diagonal with real 
positive eigenvalues Mj^^^. However, the elements of the //-matrix are, in general, complex 
quantities. Now we consider two special cases: 



(B4) 



Case-I: ji purely diagonal 



In this case, the Majorana mass matrix can be reduced to a simple block 
diagonal form which decouples the (A^i, Si) and (A^2, 'S'2) sectors: 



^ Mjvi ^ 

Mjv2 

Mn, i^n 

Mn^ /(22 J 



y2-H- r-3,c^ -H-C3 



MiVi 

Mjvi fin 

Matj 

\^ MjV2 Pi22 J 



(B5) 



18 



Then in the {Ni, Si) flavor basis, we have the 2x2 matrices 

M.= = ^\ I , (B6) 

\Mn, flu J \Mn, SiMN.e''^ 

where Si = fiu/M^. ^ 1. The Mj is diagonahzed with real and positive eigenvalues by a 
unitary transformation Uf MiUi where 

. -i cos a,e^^'/^ sina,e^^^/^ , 
Ui=\ _. _. I ' (B7) 

"sinttje '^^/^ cosOjC '^'/^ 



and the mixing angles are given by 



+ sina,~-^(^l-|) , (B8) 



up to 0{£i). The corresponding mass eigenvalues are given by 

~M^^(^1±|) (z = 1,2; J = 1,2,3,4). (B9) 

It is clear that the mass splitting within a quasi-Dirac pair is given by fiu. 

The Yukawa couplings in this diagonal mass basis are related to the couplings 
in the flavor basis as follows: 

ie-i^i/2 / £A 

hla ^ 7^ [1 + — ] Via, 



V2 \ 4 



'-2a 



V2 V'-T'^^ 



t4Q. 



v-2 O-ijf- 's^'" 

Note that in the L-conserving hmit Ei — )■ 0, we have hia = ihja within a quasi-degenerate 
pair as expected. 

Now let us calculate the CP-asymmetry for the decay of one of the quasi-Dirac 
particles, say i = 1. We have from Eq. 

/ \ 2" 



1 



Im 


'{hh% 




hi/3 


2 



- ^ V L^'^'" ''^i r Im 



i(ei-02) 



fi3 (Bll) 



assuming /13 ~ /14. Note that the j = 2 term vanishes as there is no imaginary part in 
that case. It is clear that ei vanishes as fi22 0. Similarly, one can show that €2 also 
vanishes in the limit /i22 — ^ 0, and es, €4 vanish as fiu — 0. 
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Case II: fi purely ofF-diagonal 



In this case, the Majorana mass matrix in the {Ni, Si) flavor basis reduces to 
the following block diagonal form: 



Mn^ 


Mn, /xi2 

\^ Mat, //i2 



Cl-H- C3,T'^ -H-r4 



Mn^ 

Hi2 Mn, 

Mn, 1^12 

J 



(B12) 



which, however, mixes the (1,2) sectors; in the (iVi,5'2) basis, we have the 2x2 mass 
matrix 

\ 

(B13) 



M 







with fi <^ Mn-^,Mn2- However, unlike in Case 1, we cannot diagonalize this asymmetric 
matrix by a single unitary transformation; instead, we have to apply a bi-unitary trans- 
formation of the form V^MU. We find that the following forms of U and V diagonalize 
M: 



I 



U ^ 



cos a 



sm a 



\ 



— sm ae cos ae 



V 



cos (3 sin f3 
-sin/3e^^ cos/3e^^ 



(B14) 



where the mixing angles are given by 



cos a 



sm a — 



cos /3 — 



sm / 



Ml, 













-Ml^f + ^,^Mk 


Ml, - Ml^ 


^(Mk 


-MlX + l^mk 







Mk-Mkr+i,mk 



(B15) 



The eigenvalues of M are given by 
1 



Ml, + Mk ^{Ml, - M^y + 2/.2(M2^ + Mk) 



1/2 



(B16) 



20 



up to order Note however that in the new basis, the mass matrix is still not 



/ 



diagonal and is of the form 



with eigenvalues ±Mj. This can be diagonalized 



Mi 
Mi 

with real and positive eigenvalues |Mj| by another unitary transformation: 



V 



cos ^ sm ^ 



sm ^ cos ^ 



I ■ diag (e--/^ l) 



(B17) 



We note here that in this case, unlike in case I, there is no mass splitting within the pair 
and the two quasi-Dirac RH neutrinos are exactly degenerate. This is a general result 
that the off-diagonal elements of /i do not contribute to the mass splitting within a pair; 
they just shift the eigenvalues. Hence, the splitting can be approximated by the diagonal 
elements of fi, as in Eq. ffTOj) . 

Finally, the Yukawa couplings in the mass-diagonal basis with real and positive 
eigenvalues are given in terms of the couplings in the flavor basis as follows: 



hia = 4= (cos P yia + si 



h2a 



V2 
1 



_\fi 

sm ae ?/2« i , 



cos (3 yia — sin ae y2. 



V2 

(sin /3 yi^ 



:n 

cos ae y2, 



^(sin/3 



V2 



yia + COS ae y2 



(B18) 



Note that in the L-conserving limit /i — ?■ 0, cos a, cos /3 — 1 and sin a, sin /3 — )■ 0; it is 
clear from Eqs. ( 1B18|) that in this limit, we recover the relation hia = ihja for the 
pair. 

Using Eqs. f IBlSp . it can be shown that the CP-asymmetry, Eq. ([H]), for i = 1 

becomes 

1 



sin a cos/5 ( cos^ /3 ^ — sin^ a ^ 1^27^ I /12 
\ 7 7 / 

— - |cos a sin /3 ^cos^ f3 ^ lyijl"^ — sin^ a ^ |l/27l^^ 

+ cos a cos /3 cos /3 sin/3 ^ |yi7p ~ cos a sin a ^ I2/27P 1 \ /is 



Im (e-^5:yi>,^ 
(B19) 
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which clearly vanishes in the limit yU — (as sin a, sin /3 oc /i). Similarly it can be shown 
for other channels. 

Comparing the CP-asymmetries ei in these two cases, we find that in Case I, 
the contribution within the pair vanishes and the remaining term in Eq. (IBllll which is 
proportional to /^g is highly suppressed as Mi is not quasi-degenerate with the (M3, M4) 
pair. On the other hand, in case II, the dominant contribution comes from within the 
(Mi,M2) pair which is enhanced due to large fl2- Hence, combining these results, we 
expect that in the general case with both diagonal and off-diagonal /i-entries, the dominant 
contribution to the CP-asymmetry should come from "within the pair" decay of Ni. 
We checked numerically that this is indeed the case. 
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